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ABSTRACT 


The Sensitivity of planetary waves to various initial 
conditions and thermal forcing 1S examined using a linear 
global primitive equation spectral model. 

Initial conditions are oktained by switching on a 
analytic heat source and then integrating the model equa- 
tions out to 30 days. An averaging procedure is used to 
eliminate any transient modes which remained after integra- 


tion so that the initial conditions represent steady state 


solutions. Additional integrations are performed in which 
‘errors’ are introduced into the forcing and initial 
conditions. 


Resuits of the study are examined using polar phase vs 
amplitude plots (harmonic dials) of various spherical 
harmonics. Results indicate that planetary waves are not 
sensitive to errors in the foreing or initial Conditions 
However, this lack of Sensitivity is most likely due to the 
Simplified initial conditions used in the model rather than 


to the inherent sensitivity of flanetary waves. 
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I. INTRODUCTION 

Males greateupreagetability studies (Lorenz, 1969) indi- 
cate each scale of motion has a limit toit's  predict- 
aoa lity. Small scale motions are theoretically predictable 
to an hour, synoptic scales to a few days and planetary 
scales to a few weeks. However a number of studies, Lambert 
and Merilees (1978), Baumhefner and Downey (1978), and Morse 
(1983), have shown that synoptic scale motions, rather than 
planetary scale motions, are the most accurately forecast. 
Fig. 1.£141sS taken from Daley et al (1981). It gives the 
500mb geopotential forecast error from the 6-level National 
Center for Atmospheric Research global circulation model as 
a function of wavenumber. The error for each wavenumber has 
been normalized by the natural atmospheric variance for that 
wavenumber. Note that as previously stated the synoptic 
scales(wavenumbers 4-8) are in fact the most accurately 
forecast. This result does not seem to be dependent on the 
type cf model used. Baumhefner and Downey examined a number 
of different models and found similar results. The error in 
forecasting planetary scale waves might not seem to be that 
important given that the error in the synoptic scales is 
smaller and that these are the scales of motion which 
produce most of the day to day weather changes. The impor- 
tance of an accurate planetary scale comes to light when one 
Starts to consider medium to long range forecasts (up to 10 
days). Since the planetary waves often act to steer the 
smaller synoptic scale disturbances an improved planetary 
scale forecast would presumably lead to an improved forecast 
on the synoptic scale. In addition, the planetary waves 
contain a major portion of the eddy kinetic energy in the 
atmosphere and again a better iong wave forecast would lead 


to a better synoptic scale forecast. 
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There are at least two possible reasons why planetary 


waves are not forecast as well as theory suggests. One is 


NORMALIZED RMS ERROR 





12345 678 9 10 11 i213 415 to i716 
WAVE NUMBER 





Figure 1.1 Forecast EEroOrs (fo mat, 2 
and 3 Days of Integration. 


that the model dynamics for planetary waves are inadequate. 
Given that planetary waves are guaSi-stationary in nature 
and are forced in part by differential heating and orograyhy 
it would seem reasonable to suspect that inadequacies exist 
in the dynamics of the model's forcing. Another possibility 
is that current initialization procedures are not adequate 
for planetary scales. 

It is the purpose of this study to examine the nature of 


planetary waves and to determine the sensitivity of 
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planetary waves to initial conditions and forcing. Chapter 
two of this thesis describes the model used in this study. 
The third chapter outlines the nature of the experiment. In 
the fourth chapter the results of the experiment are 
presented and in the final section these results are 


discussed and some conciusions are drawn. 


em 


II. MODEL DESCRIPTION 

The model used in this study 1S a baroclinic spectral 
transform model which was developed by Rosmond (1977). 
Simiiar formulations have been carried out by Hoskins and 
Simmons (1975) and by Bourke (1974). The model is config- 
ured to include friction and/or diabatic heating. The 
specifics of how friction and diabatic heating are included 
in the model will be discussed in the next section. The 


basic equations of the model, in Sigma coordinates, areas 
Lollews: 








>> 
96 2 LF. c+£)0-R- Fx (RT tgo4) +k FxF (2.1) 
dt a0 

> > V ‘ 
ee FoF (CHET F+ (RTD qo ae 2 (gad )+¥-F (252) 
ot a0 2 
oq . GetQg 9h 55) 
oa a I-36 
BO LS. ee ce ao 7 (2.4) 
areas aed US eS 
ae 2 RE (2.5) 
00 O 
wheres; 


C. = Vorticaty 


D - divergence 

I = Cemperaktire 

5 ~ potential temperature 

™ - surface pressure 

V = horizontal velocityi vector 

¢? - geopotential height 

Rh =) Gas) Constant 

Cp - specific heat at constant pressure 
tf —- Corioliss parameter: 


i 


- vertical coordinate (o = P/m) 


0 
¢ - vertical velocity (6 = 2 
ge in 
P- Pp 
—- R/ED 
me- frictional force 
HHewCoOntiInulecy €dudtion (Eq. 2. 3) may be rewritten by 


integrating with respect to sigma and applying the boundary 
conditions 5(0)=9(1)=0. Thus the integral of Eq. 2.3 maybe 


written as 


= -DtC (2.6) 


e7reo 
t|0 


er 


> 
where {  ) =f )do and G= V°Vq. 

The vertical velocity, ¢, may be obtained diagnostically 
Pyesubstituting Eg. 2.6 into Eq. 2.3 and integrating in the 


vertical to obtain 


6 = (D+G)o-s°(D+G) do eo 
which uses 5 (0) =0 


The first law of thermodynanics Eq. 2.4 can be written 


oT 


OT _ B.bp_,K 29 “K Od4%, Q (2.8) 
aE 7 WVIVT-0' 6a To) 4K T (SHV Va +6, 


In order to apply semi-implicit differencing it is necessary 


to divide temperature as follows: 


T = TH(o)4T'(o, X, ¥, t) (2-9) 


where T* represents an appropriately averaged temperature. 
The basic equations can be conveniently written in 


Spherical coordinates by defining the following operator: 
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2 ae ir = (2208 


Using Eqs. 2.9 and 2.10 the basic equations can be writees 


as follows: 


3 
SE a (A,B) (2< 4H 
ae = ~a(A,B)-"? (E+6+RT#q) (22 12) 
as ~a(UD" VT" )+DD"-o 65 (To) +xT(6-G-D) (2.13) 
a ee (2. 14m 
ot ad -(D+G) 
oo? = -RI | (2. 15) 
where 
A= (c+f Ute SL) (1x? ees : 
T!9q ,V1=x? 

B= rev EO: F, 

2 _U aq yaa 
25 j-x2 noe Tres 


E = (U24V?)/2(1-x?) 
U = ul(cos¥)/r. 

Vo soy (ees 7 2 

x = sin 

Y = latitude 


h = longitude 
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Pee eae loedite tne = bDasie equations used in the nodel. 
These equations are represented spectrally in the horizontal 


and finite differenced in the vertical. 


A. VERTICAL STRUCTURE 


The vertical structure of the model follows the develop- 
ment given by Arakawa and Suarez (1983). The variables are 
staggered ino so that &, D, U, V and T are carried at the 
mid-point of each layer where enc and 6 is carried at the 
top and bottom of each layer where 50 =o% . The vertical 
structure is illustrated in Fig. 2.1. The finite difference 
HOLM Of Eqs. 2.10 - 2.15 are 


Q> 
rn 


a = -a, (A,B) (2. 16) 
oPk = a (B,-A)-¥2(o, +RT*Q+E,) (2.17) 
at ae k k 
dq _ a (Zena) 
a -(G+D) : 
otk = Lg (UT! VE')T'D. + TH(G. -G-B)- (2. 19) 
art k k °k’°k 
ier.’ p | P 
aoe Peet + AU(T, - =k T. .}+0/Cp 
ho, k+1°k Pie kel k OT ge! 
: ie niee ier | (2220 
Met, -Cp(P, =P) (=Keky ok kt) a 
k 'ktl en Pe ae 
K 
Oye $ +CpT*, (E- = (2.21) 
LM 
g = 9 (G+D)-» (G.+D.)Ac. (2.22) 
et | Ke pe J _J i) 
where 
ee K 
ev = Co, 47) 
n* l+K l+te 
= aa 
Sac? k+1 2K OEE eh) 


ie 


mm ei i a i se i 


meee i i ei i al le 





wee wee eee 


mm ee i ieee eee i i i ia i i a ei ll = 





— ii iia ii i si se a ei ela eel = 


) 


LM +1 Owes Gina 


Paguirew2. 1 Vertical Structure: 
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D> 
iT 


See ey ga, (VV) /2h0 


bes k 


#{RT" (lox?) /r7} a9 -(/I-x?/r)F, 


EB = + , J . 
ce EV +d) 4 UL ee Uy 1 / 2A. 
E t ny 3G ae 
(CDRTY /r) x3 + ( =< /r)F, 
> 
G = VirVa/n 
> 
D= V-V, 
AL = See Ege? 
ce! Py PK) =1-AL 
Becowmeeeemion=- 2.22 Can be written in matrix form, such 


that the terms on the right hand sides contain all the terms 
which are to be evaluated explicitly and the left hand side 
contains those terms which are to be evaluated implicitly. 
Eg. 2.20 can be combined with an integrated finite differ- 
Pmee LOrm Of Eg. 2.15 to obtain 


6 = C TH. oe 


Meee C IS a square matrix and the other quantities are 
column vectors. The finite difference form of the surface 


pressure tendency equation(Eq. 2.14) is 


Q 


oa - (2.24) 
ae ECG) +D, Ao, 


which can be written in matrix form as 


<4 =) -N*(G+D) (2.25) 
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ner x 
where N 1s the transpose of a constant column vector. 


Sitiiiarly Eq. 2.22 Canwbeuwr recon 


OF =e ace hp) (2. 268 

ihe next to last “term wna. ee es 

ee Boe TI t8 Ay 4 OT, - 7 (2.27) 
eal - k Pee jx-1 


For the purpose of seni-implicit formulation, the tempera- 
ture is separated according to Eq. 2.9. The mean part of 


that term can be written as 


P 
4 1B, TE panne ree pote ) = M(G#D) (2. 26m 


Eqs. 2.17, 2,19, 2.23, and 2.25 ay Now be Welee on 


aD ' = 2-28 
sctv? ($'+RT#q) =a (B,-A)- = Ky ( 
a9ents = ialé (2. 30) 
aut 

ieee (2.350 
eae DoS Kop 

Fir eale a (2). 3am 
where Q = Mt KTXN: and $¢= 9 ate and Ky, and K, represent 


terns which have not been explicitly separated out. 

The semi-implicit time differencing is achieved by eval- 
uating the terms on the left hard sides of Eqs. 2.29, 2.30 
and: 2.3 lh i1nolwer eily. The remaining terms and Eq. 32 are 
evaluated explicitly using leapfrog differencing. The 


difference equations can now Le written 


n G2 rr ?: 
Dee” (c ee area (2. 33) 
=> as 2 = = « 
Dya1-At¥? (Ee Ty tRTAq, 1) t2at(Ky) 
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T _+Atob eee) ae % Daa 
me ee AtQd. 4 2At(K.) . ( ) 
+T=> _ +T> > T > (2 25) 
Ory SEND = q__,-AtN D__4-2AtN Ci, 
Now the following eguation for D canoe yetounad by Ssubsti- 
mirerng Eq. 2.34 and 2.35 into 2.33: 
2 a pe 7 2 m ras iv (2. 36) 
aD .. = B,D _ t2at(K,) -2Atv oe et Tp 
-RT#N C 
_ n 
where the matrix operator B is 
B, = At?(C QtRTN V2 +1 ae 
B. SPECTRAL FORMULATION 
Poe eqguatesons (2600) — (2.15) are represented spectrally 
in the horizontal. The variables are represented as 
follows: 
a J m m ima De oG 
mem x, o, t) = Bey Peon ie ae = (2. 33) 
J J 
ee eG 
Te See eae = 7 
where C iS some variable and aces)” See) eer 


m = zonal wavenumber 

nh = meridional index, and n-|]m] gives the number of zeros 
between the poles (-1S x $ 1) of the associated Legendre 
function 

J = truncation limit(for this study triangular truncation 
was used with J=3) 

A= (1-1) /2 nondimensional zcnal coordinate index 

(1 < 1 < 16) 
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Note that the separation is such that the coefficients ca 
are functions of time and the vertical and the spherical 
harmonic aaeae horizontal functions of space. The normali- 
zation and orthogonality properties of the i allow the 


coefficients to be obtained as follows: 


2 ee 
f i 


It = 
C = See = al 


cy axdh (2.39) 
1B gl 


_t 
4 TT 
The non-linear terms are computed using the transform method 
following Haltiner and Williams (1980). The longitudinal 
direction is treated with a Fast Fourier Transform and the 
latitudinal direction uses GauSSian Quadrature. The number 
of latitudes, N and longitudes 8 satisfy. 

N 2 33/2 +1 M2 3044 

The number of points are chose so that there will be no 
aliasing from the product terms. For this study N=4 and 
M=16. It was discovered just before this study was 
completed that number of Gaussian latitudes used to compute 
the latitudinal integrals was actually J+1 and not 3J/2 + 1. 
However, model integrations uSing J #1 latitudes are not 
gualitatively different than these using 3J/2 + 1 latitudes. 
Even the guantitative differences are small, and will occur 
mostly in the non-linear integrations. Since the cases 
examined in this study are carried out using a linear form 


of the model, one would not expect much error. 
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III. EXPERIMENT DESIGN 

Micmincentwor tings Study 1S “to Specafy a knowh analytic 
heating function on a planetary scale and then integrate the 
equations until the model atmcsphere reaches a steady or 
quasi-steady state. These steady state solutions will be 
altered and/or the heating will be changed. The equations 
will then be integrated again to determine the sensitivity 
of planetary waves to changes in the forcing and initial 
conditions 

Heating iS introduced into the model via the thermody- 
Namic equation (2.4). QO is specified to be a Legendre 


£unction of the forn: 


ol . 
(A : M¢ gy pm imA (3.1) 
Q » x, Oo) beg Pe eo 
ets specified to have the vertical structure: 
am _ p.t-o (352) 
n ea 


Only one wave form is forced at a time and the other heating 
‘amplitudes are set to zero. The constant A 1S chosen to 
give a reasonable temperature response. In most cases a 
temperature response of a few degrees is used. 
Before the heating is "turred on" the model atmosphere 

1s specified to be ina state of solid rotation 

U = cos 

ome 
such that U = 20 m/s(15.5 °%/day) at the equator. The 
Surface pressure field is chosen to balance the mean flow. 
The initial vertical temperature distribution is given by 


the standard atmosphere values consistent with the pressure 


a 


distribution given by solid rotation. A lapse rate of 0.59 
per kilcmeter is specified tc the top of the modeu 


atmosphere. 


TABLE I 
Waves Forced for Linear Cases 


Case Wave Forced 

M n 
1 1 1 
2 1 p 
3 1 3 
4 2 2 
5 2 3 
6 3 3 


In all of the experiments the model iS run in a linear 
mode. In this mode the amplitudes of all zonal wavenumbers 
except those being forced are held constant. The model is 


run in a linear mode in these experiments to allow a more 


definitive interpretation of the results. Six iinear cases 
are integrated, each case with a different forced wave 
FOEN< For example in case 1 (see Table I) zonal wavenumber 


one (m=1), meridional index one (n=1) is forced. The results 
from these cases are fairly Similar so only the longest wave 
case (case 1) is analyzed in detail. 

Two classes of cases are examined uSing the steady state 
values. In the first, the magnitude of the heating function 
is changed and the mcdel is integrated from the steady state 
values. In the second, the s&agnitude of the steady State 
values are changed and then used aS initial conditions for 
further model integrations. In the these cases the heating 
function is the same as the one originally used to bring the 


model to a steady state. 
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Due to the similarity of the six linear cases only cases 
7 S and 6 will be discussed. These three cases represent 
the three zonal wavenumbers(1, 2 and _ 3) and the two Lasic 
Meridional structures possible in the model. Cases 1 and 6 
are symmetric in their meridional structure while case 5 is 
anti-symmetric. The primary tool used to analyze these 
cases was the harmonic dial. These diagrams are polar flots 
of phase and amplitude for a given spherical harmonic. This 
plot has the advantage that the gquasi-stationary and tran- 
Silent nature of the waves are easily Seen. Fig. 4.t1isa 
good example upon which the merits of the harmonic dial can 
be explained. In this figure the phase and amplitude for 
the temperature wave for case 1 are plotted every 12 h. The 
elapsed time since tke start of the integration is indicated 
every 5 days. For an ideal case where there is a Stationary 
and a transient component one wculd expect to see a circular 
pattern which is confined to a farticular quadrant. This is 
because as the transient wave moves through the staticnary 
wave it will come inand out cf phase with the stationary 
wave. The maximum amplitude occurs when the transient and 
Staticnary waves are completely in phase and the mininun 
amplitude occurs when the waves are 180° out of phase. A 
rough estimate of the phase speed of the transient wave can 
be obtained by observing the change in phase with time as 
plotted on the diagrams. The direction of propagation can 
be determined by observing the direction of the phase change 
on the dial. A counterclockwise change in phase indicates 
eastward propagation, a clockwise change indicates westward 
propagation. In Fig. 4.1 the phase speed of the temperature 


wave at 0 = .925 can be easily eStimated. The period of the 
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wave 1s 21.5 days, thus the phase speed in degrees per day 


is 360/21.5=16.74. The direction of propagation is east- 


ee 


PHASE VS AMPLITUDE 
FOR TEMPERATURE 





Figure 4.1 Harmonic Dial for Case 1 Temperature Component. 


ward. Another important point to note here is that, given 
the uniformity of the change in phase, the phase speed of 
the transient component is constant. The harmonic dials for 
Surface pressutre, temperature and the zonal (JU) and meri- 
dional(V¥) wind compenents for the lowest model level are 
Shown in Figs. 4.2 - 4.7 . Only those components which have 
Significant amplitude are shown in these figures. In vaian 
cases the pressure and winds waves show a much more irreg- 
ular pattern than the temperature wave. These irregulari- 
ties might be due to the presence of a number of transient 
modes components with periods much less than that indicated 
by the overall pattern. Table II gives the estimated period 
(in days) and phase speed (in °/days) for the temperature, 


Surface pressure and the V and U wind component waves for 
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Figure 4.2 Harmonic Dials for Temperature 
and Surface Pressure for Case 1. 
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Figure 4.3 
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Harmonic Dials for Meridional (U) 
and Zonal(U) Wind Case 1. 
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Figure 4.4 Harmonic Dials for Temperature 
and Surface Pressure for Case 5. 
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Figure 4.5 Harmonic Dials for Meridional (JU) 
and Zonal(U) Wind Case 5. 
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Figure 4.6 


Same as Fig 4.4 for Case 6. 
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Figure 4.7 


Same as Fig 4.5 for Case 6. 
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the three cases being considered. These estimates were 
obtained from Figs. 4.1 - 4.7. Phase SpeedS were not 
calculated for the pressure and wind waves for case 5 due to 
the highly irregular nature of those dials. The phase speed 
of the temperature, pressure and wind waves are the same for 
wave component (3,3). Such is not the case for the other 
cases, though the phase speed are similar. The computed 
phase speed of the different variables for case 1 varies 
meom 19.6 °/day for P to 18.5 °%/day for JU. The other cases 
have similar phase speeds for the different components. The 
general pattern which emerges from the data is that the 
phase speed of the transient components seems to be indepen- 


dent of wavenumber and that the waves are moving at approxi- 


TABLE Il 
Period(p), and Phase Speed(c) for Transient Components 


Case ease ease 
1 5 6 
Cc Cc Cc 
23.5 i5sem 2° : 7.5 MIGET 
feet. 5 eo. ? «612.5 16.7 735 16.7 
uv 619.5 18:4 # x 725 1627 
Wee21.5 16.7 * * 7250 1627 


Mately the speed of the mean wind(15.5 °/day). To test this 
hypothesis additional integraticns were accomplished holding 
all other factors the same but increasing the speed of the 
mean wind. The results of these integrations showed an 
increase in phase speed of the component waves corresponding 


to the increase in the mean wind speed 


A. STEADY STATE CASE 


Since the purpose of this study was to see how sensitive 
the steady state solutions are to errors in the forcing and 


initial conditions it is important that the steady state 
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solutions of the model be examined to insure that they are 
indeed steady and that these sclutions are consistent with 
known theory. 

To oktain steady state solutions the basic eguations 
were integrated up to 30 days. It is obvious from Figes 
4.1-4.7 that after a rather lengthy integration the model 
atmosfhere has yet not reached a steady state and that in 
fact the transient modes identified in the last secticn are 
clearly evident thrcughout’ the period of integration. 
Though there are indications that the solutions are slowly 
converging as indicated by the begining of spiral patterns 
in Figs. 4.1-4.7 a final steady state would not Le reached 


for some time. However, by averaging over the period of the 


transient wave, it's effect can be averaged out. The 
remaining vaiues are the steady state values. the above 
procedure was performed on case 1. To insure that the part 


of the solution which remained after averaging was truly the 
steady state soluticn another integration was performed 
using the averaged values as the ainitial conditions. 
Harmonic dials were again constructed for each of the vari- 
able wave components. If the fields obtained by averaging 
were truly the steady state solutions, one would expect to 
see no phase or amplitude change with time, i.e. the 
harmonic dial would collapse to a single point. Figs. “dae 
and 4.9 contain the harmonic dials fron a 20-day integration 
initialized with the averaged values for case 1. The 
temperature wave is almost completely stationary, but the 
pressure and velocity components have small growth and move- 
ment. The fact that the pressure and wind component waves 
are not completely stationary iS another indication that 
there are transient modes in these waves which are not 
contained in the temperature wave. Conseguently, the aver- 
aging which was performed over the period of the temperature 
wave did not average out all of tae transient parts of the 


pressure and wind waves. 
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It appears from the above integrations that solutions 
obtained by averaging over the period of the transient 
component are very Ciose to steady state values. inode r 
to better understand these solutions they will be compared 
to the steady state solutions from a simplified analytic 


model. This model will be derived in the next section. 


B. THE ANALYTIC MODEL 


The predictive system of equations for planetary scale 


motion in Z=—1n (P/P) coordinates are: 


a 

eo > - uo 
D+se -2 = 0 woe 
2 (ad, eee eas - (4.2) 
me VV 5> +2ZT(2)+2T' = Q(2) 
= ZRXT (4.3) 
moot 
where 
D= V-V 
p/p 


mu o«(Z)+o'(x, y, Z) 








mee eo? LU ie Tee 
T = =5(5 +k) = a Pee Oe 
maz | <° (os a a a 
= 
feee= RI/g 
Mero ¢ lescauk 
ec =e etc et 
Meat co OZ? 


Peete! 1S the continuity equation, Ey. 4.2 is the first law 
of thermodynamics and Eq. 4.3 is the geostrophic relation. 
These eguations were obtained ty scaling the general karo- 


clinic equations (Haltiner and Williams, 1980) with: 
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LJ/a-\k <<.01 e€7-) 100 aca 2 

Where L is the horizonal scale, R is the Rossby Numberand ¢€ 
1s a rotational Froude number. Assuming a steady state 
atmosfhere (2 = 0 ) in solid rotation (U=wa cosY ) and 
droping ther'(Z) term, because it is small, Eq. 4.2 becomes: 


a9 ad! : | (4.4) 
3 


ay a7, haa = ONEZA) 


The meridional and zonal components of the geostrophic wind 


in spherical coordinates are: 


V = i do! (4 3H 
2arcos¥sin¥ 9dr 
= J (4.6) 
2aeNnsinY¥ oy 


Using Eqs. 4.5 and 4.6 the divergence in spherical coordi- 


nates can be written: 


a 2ZQcos¥ _ _ BV (4.7) 
2arsinY — fj 

Now Eg. 4.1 can be written: 

Oi) vee 1 ae ae (4.8 

= aeaaney apne =e: : 

Assuming solutions of the forn: 

g' = o(zye"™ py) (4.9) 

2 = W(Z)e"™*p™y) (4. 10) 

Q = NCZ)e*™ Py) Oe 

Eq- 4.4 and 4.8 becone: 

wimd(Z)+TW(Z) = N(Z) a 
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eno 2 Con) 
SW aeisan2s U 


Ege 4.12 and 4.13 can be combined to obtain a Single equa- 


tion for W 


me sin: ¥(wW"'-w')+Iw = N (4. 14) 


Eq- 4.14 is just a 2nd order ncnhomoyenous ordinary differ- 


ential equation for W whose general solution is of the form: 


HE ls 22 
Where ¥, Ww are the solutions to the homogeneous part of 


Ege 4.14 and Wo 


2 
timer’ ? the solution to Eg. 4.14 is 


We C_.W.+C.W as 


is the particular solution. Assuming N(Z) = 


~¢ 2 2 are 
W= C et tec en +he ae = eee! “+C,€ d e [ee 


tf 


2. = t¥1/uU-T/wa2Qsin¥ = s +q*: 
a°*= some positive number 
A = Nj/(wa*sin? ¥(a%+02)+4Pr) 


There are two possible solution cases for the homogeneous 
Bart Of Eq. 4.14. If y2 is positive the solution consists 
of an exponentially growing and decaying part. After 
applying the boundary conditions w(0)=0 and finite energy at 
Z=o, i.e-(plol2 ) = e-“]odJ]2 the solution becomes 

: pee ee 

fae = Ate ee a a) ene 
An expression for @$(z) can be obtained by substituting fq. 
4.16 into Eq. 4.13 
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2 
6(Z) = -a2Qsin®¥Ale’ “(pene A ee (4.17) 


If q* is' negative the solution will be oscillatory in Z and 
of the forn 


ae, . 2 
WZ Sacre HAC athe) spe o o (4. 18) 
where 
i = ¥P/oaen sins Vee (4.19) 


After applying the radiation boundary condition at Z=~i.e. 
w(Z) $(Z)>0 and the lower boundary condition w=0 at Z=0 Eq. 


4.18 becomes 


: ne 
W(Z) = Al (Si-u)e' ot) 2 eee oe Ay (4. 20) 
The corresponding expression for (Z)aeas 

La ee 4 
$(Z) = a’asin?WA{(Si-uye 2**H24_4¢_1 92 )e a Ly (4.21) 


Eqs. 4.16,4.17 and Eqs. 4.20,4.21 can be substituted into 
Egse- 4.9 and 4.10 to obtain exfressions for ¢and 7 Aiten 


taking the real part the q? > 0 solution are 


eaters ee eet 


2 
g' = a*Qein-vAte (non) tem “(1+a?)}sin(ma) Pm (4. 
: - a" 4.23 
Zee Ne “-e™*)cos(ma) Pm ( 
The g2 < 0 solutions are 
g' = a*gsin? vate“! ?(3sin(uZtma) (4. 24) 


ee 
tucos(uZtmA)+( lta? )e : “sin(ma)}P™ 
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. eye: 
moa aie ~ Oe ee) )cos(md) ar 


ZL/ 2 


+(e’ “sin(uZ) )sin(md) - 


If q2> Othe resfense of the atmosphere to heating is 
trapped near the surface and the vertical structure is one 
in which the waves have no tilt with height. If g2 < 0 the 
perturbation due to the heating will radiate to great 
heights and the waves will tilt westward. In both cases the 
steady state represents a balance between advection, adia- 
batic warming(cooling) and diakatic heating(cooling). The 
longitudinal and latitudinal structure, for each case, will 
be determined by the forcing specified. Theaters, 2f Waven— 
umber one is specified in the forcing term, the steady state 
solution will have a wave one fcrm. In addition, since the 
winds are geostrophic their structure will be determined by 


the structure of $'. Eq. 4.24 can be written 


o' = -A( (3e4/*sin(yz)-e"/ *yeos(uZ) deos(m)) tee 200 
Bee 2 ae 
+(5e / een ee’ “sinuz)=(1+0206 > “)sin(ma)} P™ 
which can be written in the forn | 
(G(Z)cos(mAtH(Z) )}ACY)P™CY) (4.27) 


This form allows comparison with the phase and amplitude of 

the model Legendre functions. If the model steady state 

Solutions really are steady state one would expect them to 

at least qualitatively agree with the analytic’ solutions 

derived above. The differences between the numerical nodel 

atmosphere and the analytic atmosphere are: 

1. A nondimension vertical ccordinate Z is used in the 

analytic model while sigma coordinates are used in 
humerical model, 


2. Friction is not included in the analytic model, and 
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3. F(Z) is assumed to be constant in the analytic model 
while it 1s not in the numerical model. 
T(Z) is not constant in the numerical model because f(z) is 
constant. The fact that [l(2) is not constant given TP (Zz) 
can be shown by examining the eguation for f (Z) 


aT, (4.28) 


=~ Hote eel 
ae ooh 


After converting to z2 coordinates Eq. 4.28 becomes 
2 
H? gy aT) (4 29) 


Since the lapse rate is a constant 6.5 °/Km in the model the 
second term in Eq. 4.29 1S constant but T decreases with Z 
so the first term is not constant thus [(Z) will not be 
constant inthe model atmosphere. However, since the 
decrease 1n temperature is rather small compared to the 
mean temperature [(2Z) is nearly constant so that gqualita- 
tive comparisons are still possible. in addition the 
vertical structure of the forcing in the analytic model is a 
Simple ee forme ; while in the model it ism 
“Hs 


the form O55 = C30 ° However, given that is chosen such 


that the e-folding depth of ens is the same as that for 
ex 8 then the two functions will be similar except for the 
upper levels where both functions are small and have little 
effect on the total solution. Fig. 4.10a 1s a plot ota 
numerical model forcing vertical Structumemen cme P. is 
assumed to be constant and equal to 1000mb. This 1S note 
bad assumption because the model contains no topography and 
the maximum pressure fperturbaticns are only a few millibarcs. 
This assumption also allows qualitative comparisons between 
variables in and P coordinates. The e-folding depth here 
occurs at P=/700mb. For P=700 , Z=-1in(700/1000)—-In(a7ie 


Thus, to have the same folding depth a2 must be choosen to 


4Q 


| 
| 
| 
| 
| 
2 


PRESSURE 
400.0 200.0 


600.0 


800.0 


1000.0 


0.0 0.4 0.8 1.2 1.6 2.0 
AMPLITUDE 


(a) 


PRESSURE 
800.0 600.0 400.0 200.0 0.0 


1000.0 


aa, . ley > 
0.0 0.4 0.8 lee 1.6 2.0 


AMPLITUDE 
(Db) 


ee ee ee) oe ge oe ee eee eee em 


Figure 4.10 Heating Vertical Structure 
Ya) Numerical model and (b) Analytic. 
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be™equal to 1714-ine 7s Fig. 4.170b is “a plotter gaz 

converted to P coordinates. Note that the plots are very 
Similar up to 400Omb. Figs. 4.1lia and 4.1i11tb are plotsman 
G(Z) and H(Z) vS pressure using the following estimate for 
the various parameters in Eq. 4.19: w =(20 m/s)/a, lea 

H2N2, H=10km, N2=10-4s-1, a=6.37 x 106 m and & = 7)2gum 
10-Ss-2 , Figs. 4.12 and 4.13 contain plots of the steady 
state model solutions for the meridional (V) and zonal (U) 
wind vs sigma for case 1. Since the winds are geostrophic 
the vertical structure U and VY can be compared with the 
vertical structure of $. Note the striking similarities 
between the two sets of figures. Both amplitude plots show 
a slight decrease in the lower layers and an exponential 
increase near the top of the atmosphere. Also, given Fildes 
4.14, 4.15 and 4.16, it is obvious that the model's steady 
state solutions have the longitudinal wave one form sfeci- 
fied in the case 1 forcing tern. These figures are hori- 
zontal cross-section for vorticity, meridional velocity ang 
temperature near the top and bottom of the model atmosphere. 
Contour intervals are indicated at the top of each cross- 
section. The vertical scale at the bottom of the figures is 
used to represent topography, which is not included in this 
model. It is clear that the steady model solutions are 
qualitatively consistent with the analytic solutions and 
that the method used to determine these values was a sound 
one. The question might now be asked why are the model 
solutions oscillatory in nature? Is it not possible tha 
they could be of the exponential decaying type? These ques- 
tions Can be answered by examining &q. 4.19 more closely. 
Given that 

uF 


r/wa?sin?y > 7 (4. 30) 
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Figure 4.11 Analytic Vertical Structure 
(a) Amplitude and (b) Phase. 
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Figure 4.12 Model Vertical Amplitude 
and Phase Structure for Case 1. 
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VY Velocity Structure for Case 1, 
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Eieomsolurtton will be oscillatory. In the model, wis sfpeci- 
fied so that the zonal wind is equal to 20 m/s at the 
equator so that U=wa cos = U=wa at =0 thus w=U/a. Given 
T~H2 N2, H~10 km, N22 =107-%s-!, a=6.37 x 105 m, Q =7.292 x 
10-Ss-1 one can try to solve fcr a latitude where the solu- 
tion will no longer be oscillatory. After using the above 


values Eg. 4.30 becomes 


sin¥ = 20.84 (4-31) 
which has no solution. Thus there iS no point where the 
solution would not be oscillatory. Temnomta aN None 


oscillatory solutions one would have to increase the zonal 
winc speed to almost 100 m/s. Using these steady state 
solutions it is now possible tc examine the sensitivity of 
these solutions to changes iit nC anOrGcaing and initial 


conditions. 


Sem ooNSITLIVITY EXPERIMENTS 


Using the steady state solutions from case 1 as initial 
conditions two additional experiments are carried out to 
determine the sensitivity of the planetary waves to ‘errors! 
immene LOrcing Or initial conditions. In the first experi- 
ment errors are introduced in the forcing field. aie dc real 
100% error was introduced by turning off the forcing alto- 
gether. The model is then integrated out to 20 days. Figs. 
4.17 and 4.18 contain the harmonic dials for the above inte- 
gration. Note that for the most part the once gquasi- 
Stationary solutions are now almost purely transient since 
the phase through the entire range of phase. In the previ- 
ously observed pattern the stationary plus transient nature 
of the waves was indicated by a change of phase and ampli- 


tude which occured in a particular quadrant of the dial. 
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The estimated phase speeds obtained from Figs 4.17 and 4.18 
show these transient waves to again be moving at the speed 
of the mean wind. The fact that the waves are pure tran- 
Sient now is not surprising since the heating 1s no longer 
present to £f1x the wave to a particular location. 

Tn the second experiment errors are introduced into the 


rotational part of the steady state wind by reducing those 


components by 20%. The model is again integrated out to 20 
days. Figs. 4.19 and 4.20 cecntain the harmonic dials for 
this integration. From these figures one can see that the 


only wave which is moving eastward is the temperature wave. 
All the other component waves are moving westward. The 
temperature wave has a phase speed of about 12 °/day at the 
equator which is about 25% slower than the mean wind. It is 
evident that there are a number of high frequency components 
present in the other component waves which are not present 
in the temperature wave. Since there 1s no simple Cireusias 
pattern for the wind and pressure waves it is difficult to 
obtain the period of the transient part of the solution. 
The exact phase speed of these wave is not important. The 
important thing to note here is that the waves are moving 


Slowly westward. 
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Figure 4.17 Same as Fig 4.2 fcr the Case with No Heating. 
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Figure 4.18 Same as Fig 4.3 fcr the Case with No Heating. 
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Figure 4,20 Same as_ Fig 4.3 for the Case 
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Ve. DISCUSSION ANT CONCLUSIONS 

The results of the previous section are best discussed 
in terms of normal mode theory. A number of investigators 
(Geisler and Dickinson, 1976, Daley et al, 1981, and 
Kasahara 1976) have used this method to examine the tran- 
Sient part of planetary waves. The method basicaily 
consists of solving the linearized set of equations 
describing oscillaticns of a stratified resting atmosphere 
on a spherical earth by separating the eguations into a 
latitudinal structure equation and vertical structure equa- 
je aliey ae The vertical structure equation determines the 
vertical modes of OsCiiljatron. The separation 
constants(eigenvalues) from the vertical structure equation 
can be used to find the eigenfunctions of the latitudinal 
structure equation. The latitudinal structure equation for 
a given eigenvalue is identical to a fundamental equation of 
a free surface, if the eigenvalve is replaced by the uniforn 
depth of the fluid. This eguation was first obtained by 
Laplace in his study of free oscillations of shaliow water 
over a rotating sphere. For this reason the horizontal 
structure equation is often referred to as Laplace's tidal 
equation and the eigenvalues are referred to as equivalent 
depths. The mode associated with the largest equivalent 
depth is referred to as the external or barotropic mode. 
This is because this mode has a vertical profile of hori- 
zontal divergence with the same sign throughout the atmos- 
phere. The other modes associated with smaller eguivalent 
depths are referred to as internal modes. These modes are 
oscillatory in the vertical. For a given eguivalent depth, 
Zonal wavenumber and meridional mode there exist three sets 


of eigenfrequencies and eigenfunctions. Two are eastward- 


a 


and westward-propagating gravity waves and the third isa 
westward propagating Rossby-Haurwitz type wave. The phase 
speeds of these waves are related to equivalent depth, meri- 
dional mode and zonal wavenumber. For the Rossby type wave 
the fastest phase s peed occurs with the largest equivalent 
depth and smallest meridional mode and zonal wavenumber. 
The period of this gravest, Symmetric zonal wavenumber one 
Rossby mode is almost exactly 5 days. There is observa- 
tional evidence for the existerce of this wave(Madden and 
Julian, 1972) and there is also evidence that some of the 
planetary wave error in numerical models is in part due to 
the spurious excitation of these large scale external Rossby 
modes (Daley et al, 1981). The following transient modes 
which were observed in the three different parts of this 
study are: 
1. the mode which moved at the speed of the mean wind 
in the case where the model atmosphere was being 
brought to a steady state, 
2. the mode which moved at the speed of the mean wind in 
the case where the model heating was turned off, 
3. the westward moving modes in the case where error was 
introduced into the rotational part of the wind. 
These modes may be explainable in terms of these external 
and internal Rossby modes. The mode which is moving at the 
speed of mean wind may be a slow westward moving internal 
mode which is overwhelmed by the mean wind. The westward 
moving modes in the last experiment are most likely internal 
modes as well but are probably associated with larger eyqguiv- 
alent depths than those of the earlier experiments. Since 
phase speed for a given wavenumber and meridional mode 
increases with increasing equivalent depth, these waves 
would have a greater westward phase speed and would be less 
counterbalanced by the mean wind. Given a large enough 


westward phase speed one would actually see a westward 
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propagation indicated on the harmonic dials. The low 
frequency modes in the first two parts of this study were 
most likely excited by the discontinuous nature of the 
forcing. Numerical studies by Geisler and Dickinson (1976) 
indicate that even when the forcing is gradually turned on 
that low frequency Rossby modes are excited. An analytic 
study by Clark (1972) also showed that a mode which moved at 
the speed of tne mean wind would be excited by a switch on 
of a vertical velocity or temperature disturbance. However, 
this mode was only ocne of a number of modes which would be 
excited and it was also found to decay with time. Thus, it 
is possible that the transient modes which were observed to 
move at the speed of the mear wind are in fact internal 
Rossby modes. The fact that internal Rossby modes might 
have been excited when errors were introduced into the rota- 
tional part of the wind is consistent with the results of 
Daley et al (1981) who found that errors in the rotational 
wind will cause a spurious eéxcitation of these _ modes. 
However, it must’ be pointed out that since the results of 
this study were analyzed using Sphericai harmonics rather 
than normal modes it is impossible to confirm the above 
hypothesis. The question that might now be asked 1s why 
wasn't the external Larotropic mode excited in these experi- 
ments? The studies by Geisler and Dickinson and Clark found 
that this mode would be excited. Clark found that this mode 
would in fact be the dominate mcde. The crucial difference 
between the analytic model of Clark , the numerical model of 
Geisler and Dickinson and the model in this study is that 
this model does not have a vertical wind shear while the 
aforementioned models do. This lack of vertical shear is 
the most probable reason why the barotropic mode was not 
excited. $In the absence of a mean vertical shear the 
heating specified in this model will produce low pressure at 


the surface anda high pressure aloft. The vorticities 


5 / 


associated with each will be of the same magnitude so that 
the vertically averaged vorticity vanishes and there can be 
no barotropic mode response. When a mean vertical shear is 
present there is a differential advection affect and the 
upper and lower vorticities dcnot exactly cancel, thus 
allowing the possiblity of a barotropic mode 

Given the slow phase speed of the transient modes 
excited in this study by ‘errors’ in the forcing and initial 
conditions it does not appear that planetary waves are very 
sensitive to errors. However, this apparent lack of sensi- 
tivity is most likely due to the simplified conditions used 
in the model rather than to the inherent nature of planetary 
waves. Specifically, the lack of wind shear prohibits the 
excitation of the fast external Rossby wave. The external 
mode has a fast phase speed ( 72 °/day) and if excited it 
could rapidly propagate throughout the region of integra- 
tion. 

The results obtained by heating in this study might also 
have been obtained by includirg topography along with a 
stronger surface friction tern. Future studies should 
investigate the sensitivity of planetary wave to errors in 
topography as well as heating. In addition, future studies 
will have to use a more realistic vertical wind profile. 
Due to the iack of a realistic wind profile, which prohibits 
the excitations of the external barotropic mode, any conclu- 
sions about the sensitivity of planetary waves to initial 
conditions and forcing must be delayed until more comprehen- 


Sive studies are conpleted. 
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